Spatial fluctuations in an optical parametric oscillator below threshold with an 

intracavity photonic crystal 
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We show how to control spatial quantum correlations in a multimode degenerate optical para- 
metric oscillator type I below threshold by introducing a spatially inhomogeneous medium, such 
as a photonic crystal, in the plane perpendicular to light propagation. We obtain the analytical 
expressions for all the correlations in terms of the relevant parameters of the problem and study the 
number of photons, entanglement, squeezing, and twin beams. Considering different regimes and 
configurations we show the possibility to tune the instability thresholds as well as the quantumness 
of correlations by breaking the translational invariance of the system through a photonic crystal 
modulation. 

PACS numbers: 



I. INTRODUCTION 

Quantum correlations between components of a sys- 
tem separated spatially or temporally are on the foun- 
dations of whole new technological fields like quantum 
information and computation, quantum cryptography or 
quantum-enhanced metrology . Many successful im- 
plementations have been developed in optical systems in 
continuous variables [2, 5]. An important example is the 
intensity of light far from the single photon regime. The 
optical parametric oscillator (OPO), a device in which a 
classical input (laser) beam is down-converted in a non- 
linear medium, operates in many-photon regime. In these 
devices, the nonlinearity of the quadratic crystal allows 
for light modes interactions, and, therefore, it is at the 
origin of t he g eneration of squeezing_j6|, and entan- 
glement [8l-[lQ| in fields quadratures [3, (whose spec- 
trum is continuous). Different light modes can be dis- 
tinguishable, for instance, for their polarization or fre- 
quency. Recently, special attention has been devoted 
also to spatial degrees of freedom where quantum cor- 
relations are displayed between cavity modes or parts of 
light beams [11,]. Indeed, many applications have been 
already realized with multimode light, such as optical 
switching [l2| , quantum imaffl ng jl3M l6( , metrology p/3] , 
and quantum information [l8|, [l9| . 

A very appealing possibility is to use spatial inhomo- 
geneities in broad area devices to control relevant quan- 
tum properties of the light. In related fields, like linear 
or nonlinear classical optics, it is well-known that the pe- 
riodic modulation of the refractive index leads to gaps 
in the allowed frequencies of the electromagnetic field, a 
phenomenon known as photonic band-gap [2Q|-[22|. As a 
result, these engineered media, photonic crystals (PC), 
allow to confine and guide light leading to many appli- 
cations [20I, m, . If a transverse modulation of the re- 
fractive index is considered in nonlinear cavities, this was 
predicted to inhibit modulation instabilities at similar 
wavelengths [25|, [26[ . Recently this prediction was con- 



firmed experimentally [27|, |28| and a related phenomenon 
was proposed in semiconductor microcavities [29]. Mod- 
ulation in dissipative systems leads also to the formation 
of nonlinear structures, such as different kinds of dis- 
crete cavity solitons [30|. Furthermore, the use of PC 
is also at the basis of many proposals in quantum op- 
tics. Seminal papers [111 pointed at the possibility to 
inhibit spontaneous emission in the PC band-gaps and 
this was recently observed in different experiments (31I - 
33]. As a matter of fact, the use of PC for environment 
(dissipation) engineering, stemming from the presence of 
photonic band-gaps, is the basis of intense research activ- 
ity on cavity QED. In particular, non-Mar kovian effects 
have been predicted in quantum optics with structured 
reservoirs [3J, exploring also effects on decoherence dy- 
namics and entanglement decay [35|, [36| . 

In this work we show the effect of a transversal mod- 
ulation on quantum fluctuations and correlations in a 
nonlinear device where the presence of gaps is expected 
to inhibit quantum fluctuations. We consider a photonic 
crystal optical parametric oscillator (PCOPO), that is, 
a multimode degenerate OPO with a PC in the cavity, 
as described in |26|] where the PC effect on the mod- 
ulation instability was studied and compared with the 
homogeneous case In Ref. |38l] we presented a first 
study based on the numerical analysis of quantum fluc- 
tuations using a Langeyin treatment valid both above 
and below threshold [39[. We showed that the quantum 
correlations can be tuned by means of this PC, obtaining 
noise reduction in field quadratures, robustness of squeez- 
ing in a wider angular range, and, most remarkably, an 
improvement of entanglement above threshold [38|. In 
this paper, we present analytical results valid below the 
parametric threshold and based on linear and few-modes 
approximations in good agreement with numerical simu- 
lations of the full model. We calculate intensity fluctu- 
ations and correlations as well as quadratures squeezing 
and entanglement showing the effect of the modulation. 
The homogeneous multimode OPO was shown to present 
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squeezing, entanglement and twin beams correlations be- 
tween spatial modes below threshold (4Q|-[42| and above 
threshold in presence of stable patterns or even frozen 
chaos [39I, l4p. Similar effects have been predicted in 
Kerr media [4^ and in second harmonic generation [45| 
and in the last years there have been several successful 
experimental realizations [l3l-[l6[ l46l - [49| . The effects of 
a spatial modulation here discussed for an OPO can also 
be generalized to these other nonlinear devices modified 
by the inclusion of an intracavity PC. 

The paper is organized as follows. In Sec. [Ill we present 
the model of a PCOPO using linear and few modes ap- 
proximations for the light fluctuations below threshold 
giving the output signal field in terms of the input one. 
We also introduce a set of non-linear Langevin equa- 
tions nil that we numerically simulate to test our 
approximations. In Sec. IIII) we obtain an analytical ex- 
pression for the intensities of the signal field, showing how 
the instability threshold for signal emission can be either 
raised or lowered by means of the PC spatial modulation. 
Then, in Sec. IIVI we obtain the expression for different 
quantum correlations, such as squeezing, entanglement 
and twin beams correlations. Last section is devoted to 
our conclusions. 



II. FEW MODE APPROXIMATION FOR THE 
PCOPO 

We consider a planar cavity filled with a non linear x^^^ 
medium with one of the mirrors only partially reflecting. 
The pump field at frequency 2uj is down-converted in a 
signal at frequency cj, with polarization orthogonal to 
the pump one. This constitutes an implementation of a 
type I degenerate OPO. The input beam is a plane wave 
propagating in the z-direction (the cavity axis) with am- 
plitude assumed real. Here, we consider the effect 
of the transversal inhomogeneity of the medium filling 
the cavity. This corresponds, for example, to the intro- 
duction of a planar photonic crystal (PC) with refrac- 
tive index modulation in the plane perpendicular to the 
light propagation direction. A sketch of the device is 
provided in Fig. 1 of Ref. The intracavity dynamics 
of this PCOPO can be described in terms of a continu- 
ous of boson spatial modes Ao,i(x, t) at frequencies cjq,!, 
^. These operators obey equal time commutation 



X G 
relations llll: 



i,(x,t),i](x^t) = (5,,(5(x'-x). 



(1) 



The Hamiltonian operator reads [39|, |41| 

H = Hj fd^^J2 [^1 (^) - c^V) A, 



=0,1 



lE (4 - io) + i| (ioil" - h.c.) (2) 

where the first term describes diffraction of the fields in 
the cavity, with the Laplacian in the transverse plane 



and diffraction strengths c\ = 2co. The second term 
accounts for the interaction with the external pump E 
while the nonlinear interaction between both modes is 
given by the third term, being the coupling constant g 
proportional to the second-order susceptibility x*^^^- The 
coefficient 7 is the cavity damping rate (introduced for 
convenience as a scaling). 

The main difference with respect to a generic OPO [39|, 
41] is that in a PCOPO the intracavity photonic crystal 
gives rise to a spatial modulation of the cavity detun- 
ings Ao,i(x) [38|. This constitutes a breaking of the 
translational symmetry of the system with deep conse- 
quences both in the macroscopic fields dynamics and in 
the correlations between fluctuations. In Ref. [38] numer- 
ical results about quantum effects both below and above 
threshold where reported, based on simulation of the 
quantum fields dynamics in the Q-representation. From 
the methodological point of view this description, first 
discussed in Ref. [39| , allows one to take into account the 
full nonlinear dynamics with the drawback of being not 
amenable for analytical calculations. In the following we 
introduce a simplified and approximated model that we 
use in order to obtain analytical results below the insta- 
bility threshold. 

The intracavity fields operators Aq, A\ obey the 
Heisenberg equation 



dt 



H,A, 



■ -fAj + ^2'yA 



(3) 



where the dissipative contribution characterizes such an 
open system, with y^A^" incoming quantum fluctua- 
tions [50| . Due to the cubic form of the Hamiltonian (j2j) , 
the dynamic equations for the operator moments form an 
infinite hierarchy of coupled equations, which in turn is 
unsuitable to handle analytically. 

A commonly invoked approximation based on a sys- 
tem size expansion is the linearization around a macro- 
scopic steady state leading to a dynamical evolution for 
the quantum fluctuations governed by a quadratic Hamil- 
tonian [51]. Let us identify in each field operator a ref- 
erence average value Aj and a small fluctuation around 
it, cij = Aj — Aj. The reference values Aj are the ex- 
pectation values of Aj and their evolution is obtained by 
averaging the Heisenberg equations ([3]) and by approxi- 
mating all nonlinear terms as the product of first order 
moments. This procedure leads to two classical equa- 
tions whose steady state clearly depends on the regime 
in which the PCOPO is considered. If we consider pump 
values E such that the PCOPO is below the instability 
threshold, then the signal operator expectation value is 
vanishing = 0, independently on the presence of the 
PC. The equation for the average value of the pump field 
in this regime reduces to 



dtAo 



(l + iAo(x)-zV^) Ao + £;, 



where we have introduced the scaling x' = x/y^ and 
f = 7^. This scaling is used in the remaining of this 
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article together with the scahng for the field variables 
described in |39|] . In the following we omit the primes to 
simplify notation. 

In the case of homogeneous detuning Ao{x) = Aq 
the steady state solution of this equation is immediately 
found. Then for an OPO (or a PCOPO whose modu- 
lation is only in the signal detuning) the steady state is 
homogeneous, Aq = E/{1 + iAo). 

For a non-homogeneous pump detuning the identifica- 
tion of the stationary state is generally not trivial. In the 
following, for the sake of simplicity, we consider only one 
transverse dimension and a PCOPO sinusoidal modula- 
tion such that 



Ao(x) = Ao + Mo s'm{kpx) 



(4) 



where kp is the PC wave- number. The steady state then 
satisfies, in the Fourier space. 



Mo 



(^o(^ /Cp) 



This gives rise to coupled mode equations for varying 
k. We neglect terms with \k\ > kp, a key assumption 
justified in Sec. IIIBi We then obtain for the pump field 
three non vanishing modes 



ikx 



with 



Ag(fc = 0) 



-AQ{—kp) 
E{1 



l + ikj ' 



ikl) 



1 



(5) 

(6) 
(7) 



where we assume Aq = without loss of generality. The 
steady states of the pump and (vanishing) signal fields 
are then a reference state about which the fiuctua- 
tions operators ao,i are defined. With a standard pro- 
cedure the exact Hamiltonian is approximated to one 
quadratic in these fiuctuations [42|, [5l| . Pump and sig- 
nal Heisenberg equations are actually decoupled and the 
following dynamical equation for the signal fiuctuations 
is obtained: 

Stai = -(l + iAi-z2V^)ai+Aga|; + A/277<, 

To simplify notation in the following we drop the hats 
of the operators and we denote the fiuctuation a\ as a\. 
The modulation on the signal will be similar to the pump 
one 



Ai(x) = Ai + Ml ^\x\(kpx). 



(8) 



Notice that, in general, the amplitudes of both modula- 
tion are not equal, M\ ^ Mq. Due to the presence of 
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FIG. 1: (Color online) Average far fields obtained by numeri- 
cal simulation of the full nonlinear Langevin equations. Pump 
in black solid line and signal in green dashed line, (a) cor- 
responds to the case without PC and E = 0.999, similar to 
the case with the PC affecting only the signal, represented in 
(b) for Ml = 0.5, Mo = 0, and E = 1.028. (c) shows the 
case with Mo = 0.5, Mi - 0, and E = 0.931, similar to the 
case where the PC affects both fields, represented in (d) for 
Mo = Ml =0.5 and E — 0.956. The homogeneous compo- 
nent is always p resent and harmonics appear for the signal at 
kc = v^— Ai/2 ~ 0.7 and for the pump at 2kc. The external 
pump is a 0.1% smaller than its corresponding threshold. 



the PC, the far fields fiuctuations operators ai{k) [62*] do 
not evolve independently. Different (k wave-vectors or) 
modes are dynamically coupled 

dtai{k,t) = - (1 + zAi -^i2k'^) ai{k,t) 
Ml 

+ — {ai{k + kp, t) -ai{k- kp, t)) (9) 
+ Al{nkp,t)a\{-k^nkp,t)^ J-a'^{k,t), 

n=0,±l ^ ^ 

due to the spatial modulation of both the signal detun- 
ing (Ml 7^ 0) and the pump one [through the spatial 
harmonics AQ{nkp)]. Notice that neglecting higher har- 
monics is equivalent to consider only n = 0, ±1 in Eq. ([9]). 
Still, Eq. unveils the dynamical coupling of 6 different 
modes, as discussed in App. [Aj and further approxima- 
tions are needed in order to handle this model analyti- 
cally. 



A. Numerical simulations of fully multimode and 
nonlinear dynamics 

This leads us to examine the full PCOPO model 
in order to identify the most relevant spatial modes 
in different regimes. The full dynamics can be calcu- 
lated by numerical simulation of Langevin equations (sl, 
39], obtained by mapping the full master equation for 
the PCOPO -whose system Hamiltonian is given in 
Eq. ([2])- onto an equation of motion for the Husimi quasi- 
probability distribution Q in phase space [39|, [50| . This 
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representation is then a functional of the c-number fields 
ai{x) that are used to get the expectation values of the 
operators Ai(x)[50[. In regimes where pump intensities 
are not too high, the Husimi distribution Q dynamics 
is governed by a Fokker-Planck equation, mapped in the 
following nonlinear Langevin equations for spatial depen- 
dent pump ao and signal ai fields (39| : 

dtao{x,t) = - [(l + iAo(x))-zV^]ao(x,t) + 

dtai{x,t) = - [(1 + iAi(x)) - 2zV^] + 

ao{x^t)al{x^t) -\- ^i{x^t). (10) 

with (fo additive and (fi multiplicative, phase sensitive^ 
white noises. Notice that again the effect of the PC is 
enclosed in the spatial dependence of the detunings Ao(x) 
and Ai{x). 

To study the dynamics of the system we have simu- 
lated these equations numerically (technical details about 
numerical methods are given in Ref. [HI). For the 
OPO, without PC, it is known that, for negative sig- 
nal detuning, a modulation instability develops above 
the parametric threshold [37], with critical wavenumber 
kc = Y^— Ai/2 (here we fix Ai = — 1 so that kc — 0.7). 
Below threshold, the pump is homogeneous while noisy 
precursors at the critical wave-number are observed, with 
an average intensity increasing when approaching the 
threshold [39[. In Fig. [T^ we show the average pump 
and signal far fields in the OPO, with the homogeneous 
and critical modes excited, respectively. 

We now consider the PC OPO and in particular we 
focus on the case in which the PC has a band-gap at 
the critical wave- number of the OPO 



kf) — ^kp 



(11) 



If only the signal is modulated, for Mq = and Mi = 
0.5, the PC OPO shows average far fields similar to the 
case of the OPO. On the other hand, if we introduce a 
modulation in the pump detuning, Mq ^ 0, the pump 
field develops many even harmonics 2kp^4:kp^6kp...^ as 
shown in Fig. [it and d. It is also shown that the signal 
average far field intensity remains unchanged in all cases. 

Before to proceed to the analytical evolution of the 
correlation we show here an interesting aspect of the sig- 
nal fiuctuations as obtained by numerical simulations of 
Eqs. (fTQ|) . Below threshold noisy precursors (quantum 
images) dominate the dynamics as shown in Fig. [2^: the 
"preferred" spatial periodicity corresponds to the critical 
wave- number. In the OPO the phase of this noisy pat- 
tern is not fixed and it diffuses in space in the x direction, 
being dominated by the Goldstone mode, as discussed 
in Ref. [i^. A different behavior appears when the PC 
breaks the translational symmetry leading to a modula- 
tion of the pump (Mq ^ 0, as in Fig.[2]3). Then the noisy 
pattern in the signal appears to be spatially locked and 
there are two 7r-dephased modulated modes dominating 



the pattern. The consequences in terms of correlations 
are discussed in the following sections. 



B. Few mode dynamics 

From Fig. [1] we can now see that the assumption 
ig(x) ^ ig(0) + ig(A;p)e^^^^ + ig(-fep)e-^^^^ [Eqs. (5- 
7)] is well justified as it takes into account the most rel- 
evant modes for the pump field: higher order harmonics 
at dz2A^p are much smaller than the ones at allow- 
ing neglection of the contribution of terms with \k\ > kp, 
i.e. modulations below the PC wavelength. On the other 
hand, the introduction of the modulation does not have 
any effect on the signal, whose main components are al- 
ways the modes at the critical wave- number kc = kp/2. 
Then we restrict our analysis to the modes k = ±kp/2 in 
Eq. ([9]). Within this assumption, we reduce the study of 
the PCOPO dynamics below threshold to four coupled 
operators equations that in the frequency domain read 



Lai 



where 



ai 



(ai{kc),ai{-kc),a\{-kc),a\{kc)^ 



(12) 

(13) 
(14) 

and ai{kc^uj) and a\{kc^—uj) are denoted as ai{kc) and 
a\{kc). The vector a^" is expressed in a similar manner. 
The matrix L is 

/ 1 

L = 

V 



- too 

Ml 
2 



Ml 
2 

1 — iuj 
-5* 



-S 
1 — iuj 

Ml 
2 



From Eqs. (|6]) and (O, with kp = 

^(1 



o J- 



-s 

Ml 
2 

— iuj 



\ 



(15) 



Ag(0) = S 



J 

we obtain 
z2Ai) 



2k 



l-i2Ai+MoV2' 



and 74g(±/cp) = :^f<iS^ with 



tv = 



-Mo/2 



(16) 



(17) 



1 - z2Ai 

The output fields are obtained from the input-output for- 
malism A^^^ = ^/2JAl — Ai^ [Hll , and their dynamics (in 
the frequency domain uj) is governed by 

a°^* = (2L-^ (18) 

where I is the 4x4 identity. In the following we will 
concentrate on spatial quantum effects in the signal field. 
Therefore, we can omit without ambiguity the index 1. 
To calculate different correlations in the output variables 
given the input ones, we need to obtain the inverse of 
the matrix L (the expression for this inverse is given in 
App. [B]). This formalism is at the basis of the analytical 
quantum correlation for the output fields discussed in the 
following sections. 
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FIG. 2: (Color on line) Near field evolution of the real part of 
signal ai{x,t) in OPO (a), and in PCOPO with PC affecting 
the pump (b), both 0.1% below the corresponding thresh- 
old. Space (ordinate axis) and time (abscissa) are scaled with 
diffraction length and cavity decay as mentioned in the text. 



III. INTENSITY CORRELATIONS AND THE 
PARAMETRIC THRESHOLD IN THE 
PRESENCE OF A PC 

Let us discuss the quantum effects of introducing a 
PC in an OPO below threshold, starting from the in- 
tensity (n°^*) = (a^^*'"*" a°^*) of the most intense mode 
{kc = kp/2). The analytical expressions of second order 
moments are given in Appendix O In Fig. [3^ we repre- 
sent the spectral intensity, (n°^*(/cc, a;)), whose analytical 
expression is given in Eq. ICH for different configurations 
of the PC. The effect of the PC is immediately recognized 
on the intensity spectrum: not only the intensity at each 
frequency component can be largely increased/decreased 
with respect to the case of the OPO, but also the max- 
imum of the spectrum can appear shifted (away from 
a; = 0) when the pump detuning is modulated, that is 
for Mo 7^ 0. 

The (stationary) intensity is obtained after some stan- 
dard but cumbersome calculations (Appendix [C]) and 
reads: 

(a°-*'t(^„t)a°-^(^e,t)>= (19) 
-4\Sf {4\Sf\l + - (1 + (4 + Mf)) /a, 

where the denominator a is 



a = 16\S\^\1 + - 8|5|2 (1 + (4 + 

+ (4+M2)^ 



(20) 



with S and given in Eqs. (fT6|) and (fT7j) . respectively. 

It is easy to see that this expression reduces to the one 
given in [41] in absence of the PC, when Mq = Mi = 0, 
i.e. £^^/(l — E"^). In a previous work we showed that 
the numerical results of the full model ([10]) are in a 
good agreement with this analytical expression (see Fig. 2 
of Ref. [38]), thus justifying the assumptions described 
above. 

In general, for the PCOPO, we see that both pump and 
signal modulations can modify the signal emission at the 
most intense spatial mode kc- In Fig.jSb we represent the 
intensity obtained from Eq. (p!9|) as a function of Mq and 
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FIG. 3: (Color on line) Intensity in frequency and time do- 
main for fixed values of the external pump, and different 
configurations of the PC. In (a) we show the intensity for 
E — 0.92 in frequency domain for four different configura- 
tions. In (b) we show it in the time domain for the same pump, 
and for different combinations of the parameters Mo , Mi . In 
this case, the threshold is not reached for any configuration. 
In (c) we show that, if the external pump is increased to 
E = 0.93, the intensity tends to infinity for the combinations 
of the parameters inside the white circle. 
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Ml, for a fixed value of the external pump, E = 0.92. 
The intensity can be controlled through the PC and a 
strong effect is found for small modulations of the signal 
detuning Mi and certain values of the modulation of the 
pump Mo ( Ml c:^ and Mq 0.8 in Fig. [3b). While for 
this value of the external pump E the intensity remains 
finite for every combination of Mi and Mq, this is not 
necessarily true when the pump is increased. 

Indeed, for larger pump E the average intensity in 
the PCOPO increases and, eventually, the threshold 
is reached. When approaching the parametric down- 
conversion threshold (coinciding with the spatial insta- 
bility one) our approximation would fail and this is at 
the basis of the divergences appearing in Eq. ([19]). In 
Fig. [3b we represent the intensity for an external pump 
E = 0.93 and different configurations of the PC. At the 
edges of the white circle the intensity tends rapidly to in- 
finity. Then, the presence of the PC has reduced the in- 
stability associated with the threshold to a value smaller 
than E = 0.93 for the points inside the white the circle, 
while in absence of a PC (for Mq = Mi = 0, Aq 0) the 
threshold is reached at £^ = 1. 

This is not the only possible scenario. As a matter of 
fact band-gaps are known to reduce spontaneous emis- 
sion |2Q|-[22|. Therefore, for our choice k^c = 2/cc, the 
PC is expected to cause reduction of the fluctuations 
that would lead to instability, inhibiting it. Indeed, the 
threshold can be either lowered or raised by means of the 
PC, as discussed in |38|. As the signal wavenumber kc is 
in the band-gap [25|, [26| , inhibition of the spatial insta- 
bility (lowering the threshold) is expected. On the other 
hand, the pump wavenumber is 2/cc, and therefore the 
PC introduces in the system exactly the wavenumber at 
which the instability process should occur in the pump, 
imprinting a nonlinear structure that favors the instabil- 
ity. Then a raise of the threshold can be expected. We 
find that these two competing mechanisms can increase 
or decrease the threshold depending on the relative val- 
ues of the different amplitudes of the spatial detuning 
Ao(x) and Ai(x). 

For Mo = we can easily obtain an expression for 
the increase of the threshold with Mi. In this case, the 
divergence of the intensity, when Eq. ([20]) vanishes, leads 
to the threshold expression 

£;thr = vTTTmi72F, (21) 

in accordance with the results given in [26|. This curve 
is indeed the black curve representing the threshold in 
Fig. [It. Moreover, in Fig. [Ud, the black curve represents 
the thresholds for Mq ^ and Mi = 0, thus showing 
that the threshold can be either raised or lowered with 
Mq. Finally, notice that if we increase further E eventu- 
ally reaching E = 1^ the intensity remains finite in some 
regions, in accordance with the fact that the threshold is 
increased for some configurations of the PC. 



IV. QUANTUM CORRELATIONS IN THE 
PRESENCE OF A PC 

As clear from Fig. [2] there are profound effects on the 
fluctuations of the signal field of an OPO when including 
the PC modulation. Changes in the correlations between 
the most intense modes a^^^{kp/2) and a^^^{—kp/2) are 
then expected and in the following we consider quantum 
effects such as squeezing, entanglement and twin beams 
correlations. 



A. Squeezing in the presence of PC 

Optical parametric oscillators are well-known as 
sources of spatial squeezing as mentioned in the Intro- 
duction [41, 42]. Two mode squeezing appears in quadra- 
tures of the superposition of two opposite signal modes: 

T.e^{kc,-kc) = (22) 
[a°^* {kc) + e^^a°^* {-kc)] e'^ + h.c, 

where is the quadrature angle and is a relative phase 
or superposition angle. This is equivalent to a sum of 
(position) quadratures T^o^^ {kc^ —kc) = xi + f 2, where: 

xi = a°^* {kc) e'^ + a°^*'t {kc) e''^ (23) 
X2 = a°^* {-kc) e^(^+^) + a°^*'t {-kc) 6"^^^+^) 

are the quadratures corresponding to the mode at kc and 
— /Cc, respectively. 

The variance A Ylg^p {kc-, ~kc) can be obtained in terms 
of the operator moments describing the fields out of the 
cavity as given in App. [O We remind that the first order 
moments of the signal field in all spatial modes vanish be- 
low threshold, so that /S?T^o^p is obtained from the second 
order moments of field operators. 

From the analytical expression of the variance we can 
find the squeezing for different combinations of Mq, Mi, 
and the pump E. The latter is always considered be- 
low threshold for every combination of Mq and Mi. We 
observe that the level of squeezing reached depends not 
only in the pump E but also in the values that define 
the PC. For example, for a fixed value of the pump E it 
depends on Mq and Mi, as represented in Fig. [4^. On 
the other hand, this level also depends on E and Mq 
for fixed values of Mi (see Fig. [JJd), or Mi for fixed val- 
ues of Mo (see Fig. [3t). The level of squeezing depends 
on the distance to the threshold, being enhanced when 
the system gets closer to the instability threshold. In an 
OPO, one can drive the system towards the instability 
by increasing E. A PCOPO permits one to control the 
value of the threshold through the PC, thus allowing one 
to control the squeezing by modifying Mq and Mi. In 
Fig. [4]i we represent how the variance varies with the 
quadrature angle ^, for the a value of the superposition 
angle (j) at which the maximum squeezing is achieved. 
For all four cases, we consider that the external pump is 
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FIG. 4: (Color online) Minimum value of the variance for dif- 
ferent configurations of the PC. (a) Minimum of the variance 
for an external pump below threshold, ^ = 0.92. (b) The 
variance for every for ^ at a 5% below threshold. Here, is 
fixed to the value at which the minimum variance is obtained, 
(c) Minimum value of the variance for different combinations 
of Mo and E below threshold, with when Mi = 0. (d) same 
when Ml 7^ and Mo = 0. This minimum is never smaller 
than the theoretical value in absence of PC, and it reaches 
this value at the corresponding threshold ^thr, which in (c) 
and (d) is represented by the black solid curve. 



a 5% smaller than its corresponding value at threshold 
for each configuration. Since the distance to the thresh- 
old is similar, the minimum variance reached for every 
configuration is also similar, and thus the same level of 
squeezing is obtained for different configurations. 

In order to check the validity of our approximations we 
now consider also numerical results obtained from the 
full nonlinear and multimode model. Fig. [5] shows the 
variance for the fields for all the angles and different con- 
figurations, comparing analytical and numerical results 
obtained after simulation of Eqs. ([10]), leading to intra- 
cavity correlations. The largest squeezing for an OPO 
in absence of PC is attained just at threshold between 
the critical modes while in the presence of the PC this 
minimum value is changed. The agreement between nu- 
merical and analytical calculations is very good. 

In general, we find that even considering an OPO and 
a PCOPO at a fixed distance from the respective thresh- 
olds, the variance of T^Oip is deeply modified by the pres- 
ence of the PC. While the achieved level of squeezing 
is not affected by the spatial modulation, the angles at 
which the maximum squeezing is reached are modified. 
On the other hand, we mention that the results pre- 
sented in Ref. [38| concerning the same device in the 
above threshold regime, indicate strong effects such as 
the widening in the region of angles describing maximal 
squeezed states. 



-e- 
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FIG. 5: (Color online) Variance for different configurations of 
the PC at all angles. Left column contains analytical results, 
compared with the right column of numerical ones. All the 
cases are represented for a 5% below threshold, (a) and (b) 
are for PCOPO with PC only affecting the pump (Mo = 0.5, 
Ml = 0) while (c) and (d) with PC only affecting the signal 
(Mo = 0, Ml = 0.5). 



B. Entanglement measures 

In this Section we analyze the spatial entanglement 
between opposite modes for different parameters in or- 
der to verify if the action of the refractive index modula- 
tion disturbs entanglement measures and/or creates sep- 
arable states. It is known that OPOs are useful devices 
giving spatially entangled states generation below thresh- 
old [42, 53, 54]. Let us analyze in this regime how the 
modulation introduced by the PC modifies quantum cor- 
related states at opposite points in the far field. First, we 
define the (momentum) quadrature operators pi. These 
quadratures are defined similarly to x^, Eqs. (|23j) . but 
with angles shifted as 6 — 7r/2 and (p (p -\- tt. 
Then, we introduce the following operators in order to 
establish if there is Einstein-Podolsky-Rosen (EPR) en- 
tanglement [55l] : 

. , 1. 

U \(j\Xi H X2 

(J 

V = \a\pi - -p2 
a 

where a is a real parameter. The state inseparability cri- 
terion introduced in [56] for continuous variables systems, 
establishes that any separable quantum state character- 
ized by a density operator p has a lower bound on the 
variances sum. With the notation considered here, and 
introducing the parameter a in our definitions ([22]) . the 
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FIG. 6: (Color online) Entanglement measured with inseparability and E PR- entanglement criteria for different configurations 
of the PC. (a) to (d) show the inseparability criteria taking a = 1 without PC and for the cases Mo = 0.5, Mi = 0; Mo = 0, 
Ml = 0.5; Mo = 0.5, Ml = 0.5, respectively, (e) to (h) show the EPR-entanglement criterion for the same cases. All the cases 
are represented for a 5% below threshold. White lines demarcate the regions for which EPR entanglement is found according 
to both criteria. 



inseparability criterion of Duan et al [56] reads: 

A^S^^ + A2s^+^/2,^+, > 2 (^a^ + 1) . (24) 

In Fig. [6^ to d we show this sum of variances for all 
the relevant combinations of the quadrature angle 6 and 
superposition angle (p for four different configurations. 
White lines identify the regions for which the bound 
is violated, thus corresponding to entangled states. Our 
analytical results show that this inseparability region is 
modified by a spatial modulation, even when considering 
OPO and PCOPO operating at the same distance from 
the threshold and that the quadratures and superposi- 
tion angular tolerance is slightly widened when the PC 
is introduced in the signal. 

An alternative criterion to characterize EPR entan- 
glement in continuous variable systems was proposed in 
Ref. [13. We introduce: 

A2E^^ = ((fi + Ax2)'), 

minimized by 

^ -((^1^2)) 

According to EPR criterion given in [s^ , a state is EPR- 
entangled if 

In Fig. [6^ to h we show, for four different configurations, 
the analytical calculations for this quantity. White lines 



demarcate the region for which the criterion is fulfilled 
and entanglement is found. Both measures show that 
the region where entanglement can be found is slightly 
widened if the PC is introduced in the signal. This effect 
is much relevant above threshold, as shown in (ssj . 

C. Twin beams correlations 

Finally, we present analytical results to characterize 
higher order correlations, not related in a trivial way 
to entanglement [58|. We consider twin beams corre- 
lations (go, eh characterized by a negative value of the 
normal order variance (: (n(/c) — n{—k))^ :), where, ac- 
cording to the notation introduced above, n{k) is the 
fluctuation operator {n = N — (N)) of the photon num- 
ber N{k) = {A^ {k,t)A{k,t)). As we consider Gaussian 
states, all moments can be expressed as a function of 
second order moments [50], giving: 

(: {n{k) - n{-k)f :) = {a\k)a{k))^ + {a\-k)a{-k))'^ 
+ \{a{k)a(k))\^ + \{a{-k)a{-k))\^ - 2\{a{k)a{-k))\^ 
-2\{a{k)a\-k))\^. (25) 

The non classical feature of twin beams follows from 
the negativity of the variance (: (n(/c) — n(— /c))^ :) < 
and from the corresponding singularity of the P repre- 
sentation. Notice that terms (a(/c)a(/c)), in the second 
row of the (|25]) vanish for an OPO in the absence of PC 
or when Mi = 0. This is also evident from the analytical 
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FIG. 7: (Color online) Twin beam correlations obtained from the analytical expressions, (a) Twin beam correlations for different 
configurations of the PC and E = 0.92. (b) Twin beam correlations below threshold for different values of Mi and Mo = 
and normalized with the shot noise. Black solid curve represents the value of the threshold for every value of Mi. Red curve 
represents the change from negative to positive (classical) correlation, (c) same for different values of Mo and Mi = 0. Black 
solid curve represents the value of the threshold for every value of Mo. 



correlations given in the App. O Conversely, these terms 
are large in PCOPO whenever PC is affecting the signal, 
i.e., for Ml ^ 0. Moreover, the last term {a{k)a\—k) 
vanishes, both for Mq = Mi = and for Mq = and 
Ml 7^ 0, but it can contribute to the variance negativity 
for Mo ^ 0. 

In particular, for Mq = Mi = 0, Ecl (l25l) for the output 
fields yields the simple expression [4l|, [59| 



2E^ 



(26) 



with shot noise {n{kc) -n{-kc)) = 2E^/{1-E^). There- 
fore the normalized variance is always negative and con- 
stant below threshold. This it not the case for a PCOPO 
and in Fig. [7^ we show twin beams correlations normal- 
ized by the shot noise for a fixed value of the pump and 
different PC modulations. The OPO is recovered at the 
origin, while the twin beams correlation is modified by 
the PC, becoming more classical for higher values of Mi. 
Let us consider now the case of modulation only in the 
signal detuning, that is Mq = and Mi ^ 0. We find for 
the twin beams correlations 



(: {n{kc) - ni-kc))' 



SE^{4E^ ^M^ -4) 



(4 - 4£;2 + M2)2 



(27) 



with shot noise {n{kc) - n{-kc)) = 8E'^ / {4 - 4E'^ ^ M^) . 
The normalized expression of the twin beams diverges for 
4 — 4:E^ + Ml = 0, which in turn is the expression for 
the threshold, Eq. (|2T]) (black solid curve in Fig.[7l3). We 
find that, even after normalization with the shot noise, 
the strength of these correlations is dependent on all pa- 
rameters, while in the OPO it was constant for any pump 
E. Moreover, in Fig. [Tb it is clear that the correlations 
are degraded and can become classical below threshold 
(above red line). Finally, in Fig.[7h we present normalized 
twin beams correlations when only the pump detuning is 
modulated. In this case correlations remains quantum 



for all pump strengths below threshold, even if there is a 
small reduction of this quantum effect. 

An interpretation of these results can be given consid- 
ering the modulation effects from the microscopic point 
of view. A modulation in the signal (pump) detuning 
has the main effect to couple different spatial modes as 
also clear from Eq.(j9]). In the particular case of a PC 
modulation with the periodicity (pT]) there are actually 
creation and destruction processes between photons at 
the critical tilted modes ±/Cc. Even if the nonlinearity 
gives rise to simultaneous creation (annihilation) of pho- 
tons pairs in the signal, spatial modulation of the signal 
detuning leads to photons hopping between these oppo- 
site modes. In other words there is a process of creation 
of one photon (say ai{—kc)) and destruction of one in 
the opposite mode (aj(/cc)), as can be also seen from in- 
spection of the Hamiltonian (|2]). Therefore, twin beams 
correlations, due to photon pairs emission, are present 
for small hopping rates, but when the detuning modula- 
tion of the signal is increased, the variance (|27|) becomes 
classical as the twin beams are depleted uncoherently. 
On the other hand, for Mi = and modulating the de- 
tuning of the pump field, the hopping between different 
pump harmonics is still detrimental, as it triggers differ- 
ent secondary processes besides the twin photons pairs 
generation, but has a reduced effect. 



V. CONCLUSIONS 

We have shown the effect of an intracavity PC in a 
typical device displaying quantum light spatially corre- 
lated in continuous variables, as it is the degenerate OPO 
type I. Due to the translational symmetry breaking in 
the transverse profile there are several effects also ev- 
ident in the noisy precursors locking (Fig. [2]). The PC 
modulation has a deep influence in the instability process 
and parametric threshold are both raised and lowered de- 
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pending on the configuration (Fig. [5j). As a consequence 
a modification in twin beams correlations, squeezing, sep- 
arabihty, and entanglement were expected. In order to 
analytically evaluate these quantum correlations we have 
introduced two main approximations valid below thresh- 
old, linearizing around the steady state where the signal 
field vanishes and restricting our analysis to few relevant 
harmonics (see Fig. [2]). We have considered a PC modu- 
lation (pTj) and five modes, /c = 0, ±/Cpc for the pump and 
k = ±/Cpc for the signal. Under these approximations we 
obtained good agreement between our results and the nu- 
merical solution of the full multimode nonlinear model, 
Eqs. (do]). 

In our prototype model the PC leads to a sinusoidal 
variation on the refractive index and, therefore, of the 
detunings that can affect pump, signal or both fields, de- 
pending on the values of the PC modulations Mq and Mi. 
Different configurations are described by tuning these pa- 
rameters. When the PC modulates the signal field, i.e. 
Ml 7^ 0, the instability threshold rises in agreement with 
previous works concerning single resonant cavities and 
predicting pattern inhibition [25|, [26|. We demonstrate 
that in presence of a parametric process the scenario is 
more complex and actually the threshold can also be re- 
duced, when Mo ^ 0, and the instability favored. 

Non-classical phenomena such as squeezing and entan- 
glement are very sensitive to the proximity of the insta- 
bility point. Therefore the PC has deep consequences 
and can improve correlations at a given pump with re- 
spect to the OPO, at least when threshold is lowered 
(pump modulation). Apart from this effect, all corre- 
lations have been compared in OPO and PCOPO at a 
fixed distance from threshold. Then we have found that 
squeezing (Fig. |4]), separability, and EPR-entanglement 
(Fig. [6]) are preserved both in the reached values and in 
the width of the quadrature and superposition angles re- 
gions where these phenomena appear. Finally, we have 
analytically calculated twin beams correlations (Fig|7]) 
when varying the PC modulations Mq and Mi showing 
that deeper effects are actually present in these intensity 
correlations, even at a fixed distance from the threshold, 
and that in general secondary processes degrade corre- 
lations. Besides, numerical simulations above threshold 
presented in a previous work [ssj have revealed a signifi- 
cant enhancement of squeezing and entanglement. Over- 
all, the PC allows to obtain the same quantum effects 
at a lower energy and one can enhance/avoid quantum 
properties of light, just by changing pump and/or signal 
spatial modulation in the cavity. 



Appendix A: Linear dynamics of PCOPO with any 



If we Fourier transform in the temporal variable 
Eqs. and neglect all the terms with > kp in the 



signal and pump, we obtain 



(Al) 



where 6 modes are coupled between them, namely 
ai = {ai{k),ai{k ^kp),ai{k - kp), 

a[{-k),a[{-k - kp),a\{-k + kp)^ . 

We use a compact notation and denote ai{k^uj) simply 
as ai(/c), a\{k, -u) as (2j(fe), and similarly for a^". Let us 
call r]{nkp) = —iu + (l + ^Ai + i2{k + nkp)'^) , rf'inkp) = 
-iu ^ {l - iAi - i2{k ^ nkp)'^) , S = 6'g(0), and R = 
(— Mo/2)/(l -\- ikp). Then we can write the matrix Lq as 

Ml 

hi 
2 

-6'* 

This matrix allows the dynamical description of the fiuc- 
tuations for all modes such that \k\ < kp^ with continuous 
index k. 

When considering the modulation kp = kc/2 with kc 
critical wavenumber, few intense modes are relevant and 
a reduced description can be obtained, leading to the 4x4 
matrix given in Eq. (fT5]) . 



Ml 
2 


Ml 
2 


-S 






r]{kp) 





— RS 


-s 








r]{-kp) 


rs 





-s 








Ml 
i{kp) 


Ml 
2 







Ml 








-5* 


Ml 
2 








Appendix B: Solution of the input-output equation 

The output fields dynamics is governed by the Eqs. 
(fT8|l . The inverse of L is: 

/ U{uj) V{uj) W{uj) Z{w) 

1 f -V(w) C/(w) -Z{oj) W{oj) 

D{uj) W'iio) Z'{lo) U'{uj) V'{uj) 

\-Z'{uj) W'{uj) -V'{uj) U'{uj) 



where D{uj) is the determinant 
1 1 . 



D{u:) 



-ClC2 



15^1 + 



withci = 2|S'|2(1+k2), ca = -\{Ml-A{i+Lof). Finally, 
the matrix components are 



U{uj) 

U\uj) 

V{uj) 
W{uj) 



(Cl + C2f - (1 + icj)(ci + C2) + M1C3, 

(ci + C2f - (1 + ia;)(ci + C2) - M1C3, 



-(ci +C2), 



C4+— (C1+C2), \/(cj)=C4- 

-^(Cl+C2), W'{i^) = -S\c\-rC2)^ 
-S{cif^'' + C2A^), Z\uj) = 6'*(CiA^ + C2/^*), 



where C3 = \S\ — k*) and C4 = 2c3(l — iou). 
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Appendix C: Second order moments in frequency 
and time domains 

We can use the expression of the output variables 
in terms of the input ones given in Eq. (p!8|) . to calcu- 
late different second order correlations in the frequency 
domain. For example, we can consider the correla- 
tion at different ou and the same k between a^^*'"'' and 
^out solve the expression obtained for these cor- 

relations in terms of the input variables, we should 
consider that (a^"(fc, cj) a^^'t(A;^ cj^) = ^(^ - k')S{uj - 
uo')^ while any other combination vanishes [this prop- 
erty is obtained Fourier transforming the expression 
{a'''{k,t)a'''^\k',t')) = S{k - k')5{t-t')]. For example, 
we have: 



{a-^'^\kc^uj)a-^\kc^uj')) = 
W'(uj)W{-uj') - Z\uj)Z{-uj') 
D{uj)D{-uj') 



8{^)8{uj-uj'). 



We obtain 5(uj — uj') from the Fourier transform defi- 
nition, uo' = uo due to this delta function. Notice 
that the matrix offers the term a°^*''''(/cc, —oj')^ so to ob- 
tain a°^*''''(/Cc, uj') it is necessary to change the sign of the 
terms in the matrix. Then, for the spectral intensity we 
obtain: 



(a°"*'t(fc„w)a°"*(fce,w)) = 
W'{-oj)W{uj) - Z'{-oj)Z{oj) 



(CI) 



To obtain the intensity in the time domain we should 
Fourier transform the previous expression: 



— I duj- 



To solve this integral we have to consider that the de- 
nominator shows eight poles of the form: 



±i± 



and we note that there are four in the upper part of the 
complex plane and four in the lower part. Performing 
this integral using conventional methods we obtain ([19]). 
Any other second order correlation can be obtained using 
Eq. (fT8|) as described above. For example, all the non 
vanishing terms in the expression of the variance /S^T^q^^ 
are: 



{xD = 2^{a{kc)a(kc))e^'^ 

+ 2 (at {kc) a {kc)) + [at (fe,) , {k^)] , 



{xD = (a i-kc) a i-kc)) e2^(^+^) 

+ 2 (at {-kc) a i-kc)) + [at {-k^) , at (-fe,)] , 

{xiX2) = 23? (a {kc) a {-kc)) 6^^^^+^) 

+ 23? (at {kc)a{-kc))e'^. 

Finally, the Fourier transforms of all the relevant corre- 
lations used in the definition of the variance are: 

(a°^*'t(-fc„t) a°^*(fce,t)> = 4c3C5/a, 



(a°^*(fce,t) a°^*(-A:e,t)) = (2^(-2ci)c* + csce)) /a. 



(a°^*(fee, t) a°^*(fee, t)) = - (2^(2cic; - C5C7)) /c 



(a°-*'t(fe„t)a°-*'t(-A:e,t)) = (2^*(-2cJc6 + C5C*) /a. 



(a°^*'t(fe„t) a°^*'t(fe„t)) = - {2S%24cr - csc?)) /a, 

where a and ci were defined above, C5 = (4 + ), cq - 
2 + kMi, and cr = 2k - Mi. 
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